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Abstract 

The normalization of a data cube is the ordering of the attribute values. For large multi- 
dimensional arrays where dense and sparse chunks are stored differently, proper normal- 
ization can lead to improved storage efficiency We show that it is NP-hard to compute an 
optimal normalization even for 1 x 3 chunks, although we find an exact algorithm for 1 x 2 
chunks. When dimensions are nearly statistically independent, we show that dimension- 
wise attribute frequency sorting is an optimal normalization and takes time 0(dn\og(n)) 
for data cubes of size n d . When dimensions are not independent, we propose and evaluate 
a several heuristics. The hybrid OLAP (HOLAP) storage mechanism is already 19%-30% 
more efficient than ROLAP, but normalization can improve it further by 9%— 13% for a total 
gain of 29%-44% over ROLAP. 

Key words: Data Cubes, Multidimensional Binary Arrays, MOLAP, Normalization, 
Chunking 



1 Introduction 



On-line Analytical Processing (OLAP) is a database acceleration technique used 
for deductive analysis [2]. The main objective of OLAP is to have constant-time or 
near constant-time answers for many typical queries. For example, in a database 
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Table 1 



Two tables representing the volume of sales for a given day by the experience level of the 
salesmen. Given that three cities only have experienced salesmen, some orderings (left) 
will lend themselves better to efficient storage (HOLAP) than others (right). 





<1 yrs 


1-2 yrs 


>2 yrs 




<1 yrs 


1-2 yrs 


>2 yrs 


Ottawa 






$732 


Halifax 


$43 


$54 




Toronto 






$643 


Montreal 






$450 


Montreal 






$450 


Ottawa 






$732 


Halifax 


$43 


$54 




Vancouver 


$76 


$12 




Vancouver 


$76 


$12 




Toronto 






$643 



containing salesmen's performance data, one may want to compute on-line the 
amount of sales done in Ontario for the last 10 days, including only salesmen who 
have 2 or more years of experience. Using a relational database containing sales 
information, such a computation may be expensive. Using OLAP, however, the 
computation is typically done on-line. To achieve such acceleration one can create 
a cube of data, a map from all attribute values to a given measure. In the exam- 
ple above, one could map tuples containing days, experience of the salesmen, and 
locations to the corresponding amount of sales. 

We distinguish two types of OLAP engines: Relational OLAP (ROLAP) and Mul- 
tidimensional OLAP (MOLAP). In ROLAP, the data is itself stored in a relational 
database whereas with MOLAP, a large multidimensional array is built with the 
data. In MOLAP, an important step in building a data cube is choosing a normal- 
ization, which is a mapping from attribute values to the integers used to index the 
array. One difficulty with MOLAP is that the array is often sparse. For example, 
not all tuples (day, experience, location) would match sales. Because of this sparse- 
ness, ROLAP uses far less storage. Additionally, there are compression algorithms 
to further decrease ROLAP storage requirements [3,4,5]. On the other hand, MO- 
LAP can be much faster, especially if subsets of the data cube are dense [6]. Many 
vendors such as Speedware, Hyperion, IBM, and Microsoft are thus using Hybrid 
OLAP (HOLAP), storing dense regions of the cube using MOLAP and storing the 
rest using a ROLAP approach. 

While various efficient heuristics exist to find dense sub-cubes in data cubes [7,8,9], 
the dense sub-cubes are normalization-dependent. A related problem with MOLAP 
or HOLAP is that the attribute values may not have a canonical ordering, so that 
the exact representation chosen for the cube is arbitrary. In the salesmen example, 
imagine that "location" can have the values "Ottawa," "Toronto," "Montreal," "Hal- 
ifax," and "Vancouver." How do we order these cities: by population, by latitude, by 
longitude, or alphabetically? Consider the example given in Table 1 : it is obvious 
that HOLAP performance will depend on the normalization of the data cube. A 
storage-efficient normalization may lead to better query performance. 
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One may object that normalization only applies when attribute values are not regu- 
larly sampled numbers. One argument against normalization of numerical attribute 
values is that storing an index map from these values to the actual index in the cube 
amounts to extra storage. This extra storage is not important. Indeed, consider a data 
cube with n attribute values per dimension and d dimensions: we say such a cube is 
regular or n-regular. The most naive way to store such a map is for each possible 
attribute value to store a new index as an integer from 1 to n. Assuming that indices 
are stored using logn bits, this means that n\ogn bits are required. However, array- 
based storage of a regular data cube uses &(n d ) bits. In other words, unless d = 1, 
normalization is not a noticeable burden and all dimensions can be normalized. 

Normalization may degrade performance if attribute values often used together are 
stored in physically different areas thus requiring extra 10 operations. When at- 
tribute values have hierarchies, it might even be desirable to restrict the possible 
reorderings. However, in itself, changing the normalization does not degrade the 
performance of a data cube, unlike many compression algorithms. While automati- 
cally finding the optimal normalization may be difficult when first building the data 
cube, the system can run an optimization routine after the data cube has been built, 
possibly as a background task. 



1.1 Contributions and Organization 

The contributions of this paper include a detailed look at the mathematical founda- 
tions of normalization, including notation for the remainder of the paper and future 
work on normalization of block-coded data cubes (Sections 2 and 3). In particu- 
lar, Section 3 includes a theorem showing that determining whether two data cubes 
are equivalent for the normalization problem is Graph IsOMORPHiSM-complete. 
Section 4 considers the computational complexity of normalization. If data cubes 
are stored in tiny (size-2) blocks, an exact algorithm can compute the best normal- 
ization, whereas for larger blocks, it is conjectured that the problem is NP-hard. 
As evidence, we show that the case of size-3 blocks is NP-hard. Establishing that 
even trivial cases are NP-hard helps justify use of heuristics. Moreover, the optimal 
algorithm used for tiny blocks leads us to the Iterated Matching (IM) heuristic pre- 
sented later. An important class of "slice-sorting" normalizations is investigated in 
Section 5. Using a notion of statistical independence, a major contribution (The- 
orem 18) is an easily computed approximation bound for a heuristic called "Fre- 
quency Sort," which we show to be the best choice among our heuristics when the 
cube dimensions are nearly statistically independent. Section 6 discusses additional 
heuristics that could be used when the dimensions of the cube are not sufficiently 
independent. In Section 7, experimental results compare the performance of heuris- 
tics on a variety of synthetic and "real-world" data sets. The paper concludes with 
Section 8. A glossary is provided at the end of the paper. 
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2 Block-Coded Data Cubes 



In what follows, d is the number of dimensions (or attributes) of the data cube C 
and rii, for 1 < i < d, is the number of attribute values for dimension i. Thus, C has 
size n\ x . . . x n^. To be precise, we distinguish between the cells and the indices 
of a data cube. "Cell" is a logical concept and each cell corresponds uniquely to a 
combination of values (vi,V2, . . . ,vj), with one value v ; for each attribute i. In Ta- 
ble 1, one of the 15 cells corresponds to (Montreal, 1-2 yrs). Allocated cells, such as 
(Vancouver, 1-2 yrs), store measure values, in contrast to unallocated cells such as 
(Montreal, 1-2 yrs). From now on, we shall assume that some initial normalization 
has been applied to the cube and that attribute i's values are {1,2, . . .n/}. "Index" 
is a physical concept and each J-tuple of indices specifies a storage location within 
a cube. At this location there is a cell, allocated or otherwise. (Re-) normalization 
changes neither the cells nor the indices of the cube; (Re- formalization changes 
the assignment of cells to indices. 

We use #C to denote the number of allocated cells in cube C. Furthermore, we say 
that C has density p = — x #C XWrf ■ While we can optimize storage requirements and 
speed up queries by providing approximate answers [10,11,12], we focus on exact 
methods in this paper, and so we seek an efficient storage mechanism to store all 
#C allocated cells. 

There are many ways to store data cubes using different coding for dense regions 
than for sparse ones. For example, in one paper [9] a single dense sub-cube (chunk) 
with d dimensions is found and the remainder is considered sparse. 

We follow earlier work [2,13] and store the data cube in block^i which are disjoint 
d-dimensional sub-cubes covering the entire data cube. We consider blocks of con- 
stant size mi x ... x m^; thus, there are \^~\ x . . . x \^-] blocks. For simplicity, 
we usually assume that divides n^ for all k E {1, . . . , J}. Each block can then 
be stored in an optimized way depending, for example, on its density. We consider 
only two widely used coding schemes for data cubes, corresponding respectively 
to simple ROLAP and simple MOLAR That is, either we represent the block as a 
list of tuples, one for each allocated cell in the block, or else we code the block as 
an array. For both extreme cases, a very dense or a very sparse block, MOLAP and 
ROLAP are respectively efficient. More aggressive compression is possible [14], 
but as long as we use block-based storage, normalization is a factor. 

Assuming that a data cube is stored using block encoding, we need to estimate the 
storage cost. A simplistic model is given as follows. The cost of storing a single cell 
sparsely, as a tuple containing the position of the value in the block as d attribute 
values (cost proportional to d) and the measure value itself (cost of 1), is assumed 
to be 1 + ad, where parameter a can be adjusted to account for size differences 

1 Many authors use the term "chunks" with different meanings. 
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between measure values and attribute values. Setting a small would favor sparse 
encoding (ROLAP) whereas setting a large would favor dense encoding (MOLAP). 
For example, while we might store 32-bit measure values, the number of values 
per attribute in a given block is likely less than 2 16 . This motivates setting a = 
1/2 in later experiments and the remainder of the section. Thus, densely storing a 
block with D allocated cells costs M = mi x . . . x mj, but storing it sparsely costs 
{d/2+\)D. 

It is more economical to store a block densely if (d/2+ l)D > M, that is, if 
— z^~z — > T7TTT- This block coding is least efficient when a data cube has uni- 
form density p over all blocks. In such cases, it has a sparse storage cost of d/2 + 1 
per allocated cell if p < d j\ +l or a dense storage cost of 1 /p per allocated cell 

if p > d j2 + \ ■ Given a data cube C, H(C) denotes its storage cost. We have #C < 

H(C) 

H(C) < n\ x . . . x rid- Thus, we measure the cost per allocated cell E(C) as -^r- 
with the convention that if #C = 0, then E(C) = 1. The cost per allocated cell 
is bounded by 1 and d/2 + 1: 1 < E(C) < d/2 + 1. A weakness of the model is 
that it ignores obvious storage overheads proportional to the number of blocks, 
^- x . . . x 2*-. However, as long as the number of blocks remains constant, it is rea- 
sonable to assume that the overhead is constant. Such is the case when we consider 
the same data cube under different normalizations using fixed block dimensions. 



3 Mathematical Preliminaries 



Now that we have defined a simple HOLAP model, we review two of the most im- 
portant concepts in this paper: slices and normalizations. Whereas a slice amounts 
to fixing one of the attributes, a normalization can be viewed as a tuple of permu- 
tations. 



3.1 Slices 



Consider an n-regular J-dimensional cube C and let C,-, i d denote the cell stored 

at indices {i\ , . . . , e { 1 , . . . , n) . Thus, C has size n d . The slice C{, of C, for index 
v of dimension j (1 < j <d and l<v<n)isa<i— 1- dimensional cube formed as 

Cti u ...jj_ u i j+U ...j d = Q u ..^ij_ uV j j+u ...j d (See Figure 1). 

For the normalization task, we simply need know which indices contain allocated 
cells. Hence we often view a slice as a d — 1 - dimensional Boolean array C{. 
For example, in Figure 1, we might write (linearly) C\ = [0, 1,0,5,9,2,4,0,0] and 
C\ = [0, 1,0, 1, 1, 1, 1,0,0], if we represent non-allocated cells by zeros. Let #C J V 
denote the number of allocated cells in slice C J V . 
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Fig. 1. A 3 x 3 x 3 cube C with the slice C\ shaded. 
3.2 Normalizations and Permutations 



Given a list of n items, there are n\ distinct possible permutations noted T n (the 
Symmetry Group). If y e T n permutes i to j, we write y(z') = j. The identity per- 
mutation is denoted l. In contrast to previous work on database compression (e.g., 
[4]), with our HOLAP model there is no performance advantage from permuting 
the order of the dimensions themselves. (Blocking treats all dimensions symmet- 
rically.) Instead, we focus on normalizations, which affect the order of each at- 
tribute's values. A normalization 71 of a data cube C is a J-tuple (yi,...,y^) of 
permutations where y e T n for i = 1, . . . ,d, and the normalized data cube n(C) 
is n(C)i h ..j d = Cftfa),...^) for all (h,. . . ,i d ) e {1,. . . ,n} d . Recall that permuta- 
tions, and thus normalizations, are not commutative. However, normalizations are 
always invertible, and there are (nl) d normalizations for an n-regular data cube. 
The identity normalization is denoted / = (i, . . . ,i); whether / denotes the identity 
normalization or the identity matrix will be clear from the context. Similarly may 
denote the zero matrix. 

Given a data cube C, we define its corresponding allocation cube A as a cube with 
the same dimensions, containing O's and l's depending on whether or not the cell 
is allocated. Two data cubes C and C ', and their corresponding allocation cubes A 
and A', are equivalent (C ~ C') if there is a normalization n such that n(A) = A'. 

The cardinality of an equivalence class is the number of distinct data cubes C in this 
class. The maximum cardinality is (nl) d and there are such equivalence classes: 
consider the equivalence class generated by a "triangular" data cube Ci u ...j d = 1 if 
h < h < ■ ■ ■ < id and otherwise. Indeed, suppose that C yi{il ),...,y d (i d ) = c i l (h),-,i d (i d ) 
for all i\,...,i d , then yi(ii) <fii}i) < ■■■ <Jd(id) if and only ify^z'i) < i 2 {i 2 ) < 
• • • < y'd(id) which implies that y = yj for i e {1, ... ,d}. To see this, consider the 
2-d case where Yi(i'i) < y2{h) if and only if ^(z'i) < y'jih)- In this case the result 
follows from the following technical proposition. For more than two dimensions, 
the proposition can be applied to any pair of dimensions. 
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Proposition 1 Consider any 7i,72,7i,72 £ T„ satisfying yi(i) < 72(7) 7^(1) < 
72(7) / or a ^ 1 — 7 — n - Yi = Yi an d 72 = 72- 



PROOF. Fix z, then let k be the number of j values such that 72 (7) > 7i (z) . We have 
that 71 (/) = n — k+l because it is the only element of { 1 , . . . , n} having exactly k 
values larger or equal to it. Because 71 (z) < 72(7) ^ Yj (i) < YO')' i\{i) =n — k+\ 
and hence Y = 71 . Similarly, fix j and count z values to prove that Y2 = 72. □ 



However, there are singleton equivalence classes, since some cubes are invariant 



under normalization: consider a null data cube C, 



d 



for all (i\,...,i d ) e 



To count the cardinality of a class of data cubes, it suffices to know how many 
slices Cv of data cube C are identical, so that we can take into account the invari- 
ance under permutations. Considering all n slices in dimension r, we can count the 
number of distinct slices d r and number of copies n r \,. . .,n r j r of each. Then, the 



number of distinct permutations in dimension r is 



and the cardinality 



n r ,\\x...,xn rtdr \ 

Of a given equivalence class is ]lr=i ( n rJ !x ."'x^ i ) • ^ or exam pl e > the equivalence 

1 

has a cardinality of 2, despite having 4 possible 



class generated by C = 
normalizations. 
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To study the computational complexity of determining cube similarity, we define 
two decision problems. The problem Cube Similarity has C and C' as input 
and asks whether C ~ C'. Problem Cube Similarity (2-d) restricts C and C' 
to two-dimensional cubes. Intuitively, Cube Similarity asks whether two data 
cubes offer the same problem from a normalization-efficiency viewpoint. The next 
theorem concerns the computational complexity of Cube Similarity (2-d), but 
we need the following lemma first. Recall that (71,72) is the normalization with the 
permutation 71 along dimension 1 and 72 along dimension 2 whereas (7i,72)(/) is 
the renormalized cube. 

Lemma 2 Consider the nxn matrix I' = (71 , 72) (/). Then I' = I 71 = 72. 

We can now state Theorem 3, which shows that determining cube similarity is 
Graph IsOMORPHiSM-complete [15]. A problem II belongs to this complexity 
class when both 

• n has a polynomial-time reduction to Graph Isomorphism, and 

• Graph Isomorphism has a polynomial-time reduction to II. 

Graph IsOMORPHiSM-complete problems are unlikely to be NP-complete [16], 
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yet there is no known polynomial-time algorithm for any problem in the class. This 
complexity class has been extensively studied. 



Theorem 3 Cube Similarity (2-d) is Graph IsoMORPHiSM-comptoe. 



PROOF. It is enough to consider two-dimensional allocation cubes as 0- 1 matri- 
ces. The connection to graphs comes via adjacency matrices. 

To show that Cube Similarity (2-d) is graph isoMORPHiSM-complete, we 
show two polynomial-time many-to-one reductions: the first transforms an instance 
of Graph Isomorphism to an instance of Cube Similarity (2-d). 

The second reduction transforms an instance of Cube Similarity (2-d) to an 
instance of GRAPH ISOMORPHISM. 

The graph-isomorphism problem is equivalent to a re -normalization problem of 
the adjacency matrices. Indeed, consider two graphs G\ and G2 and their adjacency 
matrices M\ and M2. The two graphs are isomorphic if and only if there is a permu- 
tation y so that (y,y)(Mi) = M2. We can assume without loss of generality that all 
rows and columns of the adjacency matrices have at least one non-zero value, since 
we can count and remove disconnected vertices in time proportional to the size of 
the graph. 

We have to show that the problem of deciding whether y satisfies (y,y)(Afi) = M2 
can be rewritten as a data cube equivalence problem. It turns out to be possible by 
extending the matrices M\ and M%. Let / be the identity matrix, and consider two 

Ail I 

allocation cubes (matrices) A\ and A2 and their extensions A\ = 



7/0 
/ 00 



and 



A 2 = 



A 2 I I 
I 10 
I 00 



Consider a normalization n satisfying 71 (A 1) = A 2 for matrices A i,A2 having at least 
one non-zero value for each column and each row. We claim that such a n must be 
of the form n = (Yi, Y2) where yi = Y2. By the number of non-zero values in each 
row and column, we see that rows cannot be permuted across the three blocks of 
rows because the first one has at least 3 allocated values, the second one exactly 2 
and the last one exactly 1. The same reasoning applies to columns. In other words, 
if x e [7, jn] , then y(x) e [j, jn] for j = 1 , 2, 3 and i = 1 , 2. 

Let Yi j 7 denote the permutation y restricted to block j where j = 1,2,3. Define 



8 



yj = Ji\j ~ J n f° r j = 1,2,3 and i = 1,2. By Lemma 2, each sub-block consisting 
of an identity leads to an equality between two permutations. From the two identity 
matrices in the top sub-blocks, for example, we have that y} = y\ and yj = y^. From 
the middle sub-blocks, we have y\ — y\ and y\ — y\, and from the bottom sub- 
blocks, we have y^ = y\. From this, we can deduce that y} = y\ — y\ = y\ so that 
y} = y\ and similarly, y 2 = y\ and y\ = y\ so that yi = y 2 . 

So, if we set Ai = M\ and A 2 = M 2 , we have that G\ and G 2 are isomorphic if and 
only if A\ is similar to A 2 . This completes the proof that if the extended adjacency 
matrices are seen to be equivalent as allocation cubes, then the graphs are isomor- 
phic. Therefore, we have shown a polynomial-time transformation from GRAPH 
Isomorphism to Cube Similarity (2-d). 

Next, we show a polynomial-time transformation from Cube Similarity (2-d) 
to Graph Isomorphism. We reduce Cube Similarity (2-d) to Directed 
Graph Isomorphism, whichis in turn reducible to Graph Isomorphism [17,18]. 



Given two 0-1 matrices M\ and M 2 , we want to decide whether we can find (yi, Y2) 
such that (Yi, Y2)(Mi) = M 2 . We can assume that M\ and M2 are square matrices 
and if not, pad with as many rows or columns filled with zeros as needed. We 
want a reduction from this problem to Directed Graph Isomorphism. Con- 



sider the following matrices: M\ 



Afi 


and M 2 = 


0M 2 











. Both Mi and M 2 



can be considered as the adjacency matrices of directed graphs G\ and G 2 . Suppose 
that the graphs are found to be isomorphic, then there is a permutation y such that 
(y,y)(Mi) = M 2 . We can assume without loss of generality that y does not permute 
rows or columns having only zeros across halves of the adjacency matrices. On 
the other hand, rows containing non-zero components cannot be permuted across 
halves. Thus, we can decompose y into two disjoint permutations y 1 and y 2 and 
hence (y 1 , y 2 ) (Mi ) = M 2 , which implies Mi ~ M 2 . On the other hand, if Mi ~ M 2 , 
then there is (y 1 ,'/ 2 ) such that (y 1 ,y 2 )(Mi) = M 2 and we can choose y as the direct 
sum of y 1 and y 2 . Therefore, we have found a reduction from Cube Similarity 
(2-d) to Directed Graph Isomorphism and, by transitivity, to Graph Iso- 
morphism. 



Thus, Graph Isomorphism and Cube Similarity (2-d) are mutually reducible 
and hence Cube Similarity (2-d) is Graph IsoMORPHiSM-complete. □ 



Remark 4 If similarity between two n x n cubes can be decided in time cn k for 
some positive integers c and k > 2, then graph isomorphism can be decided in 
0(n k ) time. 

Since Graph Isomorphism has been reduced to a special case of Cube Simi- 
larity, then the general problem is at least as difficult as Graph Isomorphism. 
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Yet we have seen no reason to believe the general problem is harder (for instance, 
NP-complete). We suspect that a stronger result may be possible; establishing (or 
disproving) the following conjecture is left as an open problem. 

Conjecture 5 The general Cube Similarity problem is also Graph Isomor- 
PUlSM-complete. 



4 Computational Complexity of Optimal Normalization 

It appears that it is computationally intractable to find a "best" normalization n (i.e., 
71 minimizes cost per allocated cell 2s (71(C))) given a cube C and given the blocks' 
dimensions. Yet, when suitable restrictions are imposed, a best normalization can 
be computed (or approximated) in polynomial time. This section focuses on the 
effect of block size on intractability. 

4.1 Tractable Special Cases 

Our problem can be solved in polynomial time, if severe restrictions are placed 
on the number of dimensions or on block size. For instance, it is trivial to find a 
best normalization in 1-d. Another trivial case arises when blocks are of size 1, 
since then normalization does not affect storage cost. Thus, any normalization is 
a "best normalization." The situation is more interesting for blocks of size 2; i.e., 
which have m ; - = 2 for some 1 < i < d and mj = 1 for 1 < j < d with i ^ j. A best 
normalization can be found in polynomial time, based on weighted-matching [19] 
techniques described next. 

4.1.1 Using Weighted Matching 

Given a weighted undirected graph, the weighted matching problem asks for an 
edge subset of maximum or minimum total weight, such that no two edges share 
an endpoint. If the graph is complete, has an even number of vertices, and has only 
positive edge weights, then the maximum matching effectively pairs up vertices. 

For our problem, normalization's effect on dimension k, for some 1 < k < d, corre- 
sponds to rearranging the order of the n^ slices C\, where 1 < v < n^. In our case, 
we are using a block size of 2 for dimension k. Therefore, once we have chosen two 
slices C\ and C^, to be the first pair of slices, we will have formed the first layer of 
blocks and have stored all allocated cells belonging to these two slices. The total 
storage cost of the cube is thus a sum, over all pairs of slices, of the pairing-cost 
of the two slices composing the pair. The order in which pairs are chosen is ir- 
relevant: only the actual matching of slices into pairs matters. Consider Boolean 



10 



vectors b = C\ and b' = Cy. If both b ; and b',- are true, then the i th block in the pair 
is completely full and costs 2 to store. Similarly, if exactly one of b, and b', is true, 
then the block is half-full. Under our model, a half-full block also costs 2, but an 
empty block costs 0. Thus, given any two slices, we can compute the cost of pairing 
them by summing the storage costs of all these blocks. If we identify each slice with 
a vertex of a complete weighted graph, it is easy to form an instance of weighted 
matching. (See Figure 2 for an example.) Fortunately, cubic-time algorithms exist 
for weighted matching [20], and is often small enough that cubic running time 
is not excessive. Unfortunately, calculating the n^{n^ — l)/2 edge weights is ex- 
pensive; each involves two large Boolean vectors with ^]lf=i w ' elements, for a 

total edge-calculation time of © (%nf=i n i) ■ Fortunately, this can be improved for 
sparse cubes. 
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and r2 



10 1 



In the 2-d case, given any two rows, for example r\ 
then we can compute the total allocation cost of grouping the two together as 
2(#n +#r2 — benefit) where benefit is the number of positions (in this case 1) where 
both r\ and have allocated cells. (This benefit records that one of the two allo- 
cated values could be stored "for free," were slices r\ and paired.) 



According to this formula, the cost of putting r\ and r2 together is thus 2(2 + 2 — 
1) = 6. Using this formula, we can improve edge-calculation time when the cube 
is sparse. To do so, for each of the % slices C\, represent each allocated value by 
a J-tuple i - 2, . . . , ijfc_i , ijfc+i, • • • , id, ik) gi vm g its coordinates within the slice and 
labeling it with the number of the slice to which it belongs. Then sort these #C 
tuples lexicographically, in 0(#Clog#C) time. For example, consider the following 
cube, where the rows have been labeled from tq to r$ ( r ; corresponds to Cj): 



r 

n l l o l 
n 1 

r 3 1 1 
r 4 1 
r 5 1 1 



We represent the allocated cells as {(0,n), (l,n), (3,n), (0,r2), (1,^3), (2, 7*3), 
(l,f"4), (0, r 5 ), and (3,r 5 )}. We can then sort these to get (0, n), (0,r2), (0,r 5 ), 
(l,n), (1,7*3), (l 5 r 4)> (2,7*3), (3,n), (3, r 5 ). This groups together allocated cells 
with corresponding locations but in different slices. For example, two groups are 
((0,n), (0,r2), (0, r$)) and ((l,ri), (1, r^), (1,^)). Initialize the benefit value asso- 
ciated to each edge to zero, and next process each group. Let g denote the number 
of tuples in the current group, and in 0(g 2 ) time examine all (|) pairs of slices 
(51,52) in the group, and increment (by 1) the benefit of the graph edge (51,52). In 
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Fig. 2. Mapping a normalization problem to a weighted matching problem on graphs. Rows 
are labeled and we try to reorder them, given block dimensions 2x1 (where 2 is the vertical 
dimension). In this example, optimal solutions include r ,r l ,r 2 ,r 3 and r 2 ,r 3 ,ri,r . 

our example, we would process the group ((0,n), (0,r2), (0,^)) and increment 
the benefits of edges (n,r 2 ), (r 2 ,r 5 ), and (r u r 5 ). For group ((l,ri), (l,r 3 ), (l,r 4 )), 
we would increase the benefits of edges (r^r^,), {r\,r^), and (^3^4). Once all #C 
sorted tuples have been processed, the eventual weight assigned to edge (v,w) is 
2(#Q + #C^, — benefit(v, w) ) . In our example, we have that edge (ri , r 2 ) has a ben- 
efit of 1, and so a weight of 2{#r\ + #r 2 — benefit) = 2(3 + 1 — 1) = 6. 

A crude estimate of the running time to process the groups would be that each 
group is O(nfc) in size, and there are 0(#C) groups, for a time of 0(#Cn|). It can 
be shown that time is maximized when the #C values are distributed into ftC/n^ 
groups of size n^, leading to a time bound of 0(#Cn^) for group processing, and an 
overall edge-calculation time of #C{n^ + log#C). 

( k—l—i 

Theorem 6 The best normalization for blocks of size lx...x 1x2x1. ..xl can 
be computed in 0{n^ x (n\ x n 2 x . . . x nj) + ny time. 

The improved edge-weight calculation (for sparse cubes) leads to the following. 

j k—l—i 

Corollary 7 The best normalization for blocks of size 1 x ...x 1 x2x 1... x 1 can 
be computed in 0(#C(n/ ( + log#C) + ny time. 

For more general block shapes, this algorithm is no longer optimal but nevertheless 
provides a basis for sensible heuristics. 

4.2 An NP-hard Case 

In contrast to the 1 x 2-block situation, we next show that it is NP-hard to find 
the best normalization for 1 x 3 blocks. The associated decision problem asks 
whether any normalization can store a given cube within a given storage bound, 
assuming 1x3 blocks. We return to the general cost model from Section 2 but 
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choose a = 1/4, as this results in an especially simple situation where a block with 
three allocated cells (D = 3) stores each of them at a cost of 1, whereas a block 
with fewer than three allocated cells stores each allocated cell at a cost of 3/2. 

The proof involves a reduction from the NP-complete problem Exact 3-Cover (X3C), 
a problem which gives a set S and a set T of three-element subsets of S. The ques- 
tion, for X3C, is whether there isaf'CT such that each s E S occurs in exactly 
one member of T' [17]. 

We sketch the reduction next. Given an instance of X3C, form an instance of our 
problem by making a | T | x |5| cube. For s E S and T ET , the cube has an allocated 
cell corresponding to (T,s) if and only if s E T . Thus, the cube has 3\T | cells that 
need to be stored. The storage cost cannot be lower than 9 ^~^ and this bound 
can be met if and only if the answer to the instance of X3C is "yes." Indeed, a 
normalization for 1 x 3 blocks can be viewed as simply grouping the values of an 
attribute into triples. Suppose the storage bound is achieved, then at least |5| cells 
would have to be stored in full blocks. Consider some full block and note there are 
only 3 allocated cells in each row, so all 3 of them must be chosen (because blocks 
are 1 x 3). But the three allocated cells in a row can be mapped to a T E T . Choose 
it for T'. None of these 3 cells' columns intersect any other full blocks, because 
that would imply some other row had exactly the same allocation pattern and hence 
represents the same T, which it cannot. So we see that each s E S (column) must 
intersect exactly one full block, showing that T ' is the cover we seek. 

Conversely, suppose 1' is a cover for X3C. Order the elements in l' arbitrarily 
as To, 7\, . . . , 7|5|/3 and use any normalization that puts first (in arbitrary order) the 
three s E Tq, then next puts the three s E T\, and so forth. The three allocated cells 
for each 7] will be together in a (full) block, giving us at least the required "space 
savings" of = \S\. 

Theorem 8 It is NP-hard to find the best normalization when 1x3 blocks are used. 

We conjecture that it is NP-hard to find the best normalization whenever the block 
size is fixed at any size larger than 2. A related 2-d problem that is NP-hard was 
discussed by Kaser [21]. Rather than specify the block dimensions, this problem 
allows the solution to specify how to divide each dimension into two ranges, thus 
making four blocks in total (of possibly different shape) . 



5 Slice-Sorting Normalization for Quasi-Independent Attributes 

In practice, whether or not a given cell is allocated may depend on the correspond- 
ing attribute values independently of each other. For example, if a store is closed 
on Saturdays almost all year, a slice corresponding to "weekday=Saturday" will be 



13 



sparse irrespective of the other attributes. In such cases, it is sufficient to normalize 
the data cube using only an attribute-wise approach. Moreover, as we shall see, one 
can easily compute the degree of independence of the attributes and thus decide 
whether or not potentially more expensive algorithms need to be used. 

We begin by examining one of the simplest classes of normalization algorithms, 
and we will assume n-regular data cubes for n > 3. We say that a sequence of 
values xi, . . . ,x M is sorted in increasing (respectively, decreasing) order if X[ < Xj+i 
(respectively, Xi > for i G { 1 , . . . , n — 1 } . 

Recall that C J V is the Boolean array indicating whether a cell is allocated or not in 
slice C J V . 

Algorithm 1 (Slice-Sorting Normalization) Given an n-regular data cube C, then 
slices have n d ~ l cells. Given a fixed function g : {true, false} 11 — > R, then for each 
attribute j, we compute the sequence fi =g(C J v ) for all attribute values v= l,...,n. 
Let y 7 be a permutation such that y J (f J ) is sorted either in increasing or decreasing 
order, then a slice-sorting normalization is (y 1 , . . . , n f). 

Algorithm 1 has time complexity 0(dn d + dnlogn). We can precompute the aggre- 
gated values fv and speed up normalization to 0(dn\og(n)). It does not produce a 
unique solution given a function g because there could be many different valid ways 
to sort. A normalization 03 = (y 1 , . . . , f 1 ) is a solution to the slice-sorting problem if 
it provides a valid sort for the slice- sorting problem stated by Algorithm 1 . Given 
a data cube C, denote the set of all solutions to the slice-sorting problem by Sc,g- 
Two functions g\ and g2 are equivalent with respect to the slice-sorting problem if 
Sc.gi = ~Sc.g 2 f° r ai l cubes C and we write g\ ~ g2 . We can characterize such equiv- 
alence classes using monotone functions. Recall that a function h : R — >■ R is strictly 
monotone nondecreasing (respectively, nonincreasing) if x < y implies h(x) < h(y) 
(respectively, h(x) >h(y)). 

An alternative definition is that h is monotone if, whenever x\, . . . ,x n is a sorted list, 
then so is h(x\), . . . , h(x n ). This second definition can be used to prove the existence 
of a monotone function as the next proposition shows. 

Proposition 9 For a fixed integer n>3 and two functions COi , (02 : D — > R where D 
is a set with an order relation, if for all sequences x\, . . . ,x n G (D, C0i(xi), . . . ,t0i(x M ) 
is sorted if and only if (02(^1), • • • , 0>i(x n ) is sorted, then there is a monotone func- 
tion h : R — > R such that (Oi = h o CO2. 

PROOF. The proof is constructive. Define h over the image of (02 by the formula 

h(&2(x)) = t0i(x). 

To prove that h is well defined, we have to show that whenever (02(^1) = ©2(^2) 
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then Cfli (jq ) = (Oi (^2) • Suppose that this is not the case, and without loss of general- 
ity, let (Oi (x\ ) < (Oi (x2). Then there is X3 G "D such that a>i (x\ ) < a>\ (^3) < C0i (^2) 
or Cfli(x3) < tOi(xi) or (01(^2) < G>i(*3)- In all three cases, because of the equal- 
ity between (02(^1) and 0)2(^2), any ordering of 0)2(^1), 0)2(^2), ©2(^3) is sorted 
whereas there is always one non-sorted sequence using a>i . There is a contradic- 
tion, proving that h is well defined. 

For any sequence x\,X2, *3 such that (02(^1) < ©2(^2) < ©2(^3), then we must either 
have Cfli(xi) < C0i(x2) < 0)1(^3) or (£>i(x\) > Q>i(x2) > 0)1(^3) by the conditions of 
the proposition. In other words, for x < y < z, we either have h(x) < h(y) < h(z) or 
h(x) >h(y) > h(z) thus showing that h must be monotone. □ 

Proposition 10 Given two functions g\,g2 '■ {truejalse] s — > IR, we have that 

Sc, gi = Sc, 82 

for all data cubes C if and only if there exist a monotone function ^ : IR — > IR such 
that g\ =hog 2 . 

PROOF. Assume there is h such that g\ = h o g 2 , and consider 03 = (y 1 , . . . , y**) G 
Sc,gi for any data cube C, then yi(gi (Cl)) is sorted over index v G {1, . . . ,n} for 
all attributes j = 1, . . . ,n by definition of Sc. gl ■ Then {h(g\(Cl))) must also be 
sorted over v for all j, since monotone functions preserve sorting. Thus 03 G Sc, g2 - 

One the other hand, if Sc,gi = Sc, S2 f° r ai l data cubes C, then h exists by Proposi- 
tion 9. □ 



A slice-sorting algorithm is stable if the normalization of a normalized cube can be 
chosen to be the identity, that is if 05 G Sc, s then / G State), g for all C. The algorithm 
is strongly stable if for any normalization 03, S®(c).g 03 = 5c.g for all C. Strong 
stability means that the resulting normalization does not depend on the initial nor- 
malization. This is a desirable property because data cubes are often normalized 
arbitrarily at construction time. Notice that strong stability implies stability: choose 
03 G Sc, g - Then there must exist £ G S®(c),g such that C, o 03 = 03 which implies that 
£ is the identity. 

Proposition 11 Stability implies strong stability for slice-sorting algorithms and 
so, strong stability stability. 

PROOF. Consider a slice-sorting algorithm, based on g, that is stable. Then by 
definition 

aescj^ieSanQj (l) 
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for all C. Observe that the converse is true as well, that is, 

I e Sm(C),g 03 e s c,g- (2) 

Hence we have that G5i o 03 G Sc.g implies that / G 5 roi ( ro (c)),g by Equation 1 and 
so, by Equation 2, 03i G 5 ra (c),g- Note that given any 03, all elements of 5c,g can be 
written as 03i o03 because permutations are invertible. Hence, given 03i o 03 G Sc,g 
we have CDi G S^ C ),g and so ^C,g C S®( C ),g ° ®- 

On the other hand, given 03 1 o 03 G S®(c),g ° ®> we nave that 03 1 G S®(c),g by cancel- 
lation, hence / G ^(osfc)),,? by Equation 1, and then 03i o 03 G Sc.g by Equation 2. 
Therefore, S®( C ),g ° 05 c ^c,g- 1=1 

Define x : {truejalse} s — > R as the number of frw values in the argument. In effect, 
x counts the number of allocated cells: x(Cy) = #C J V for any slice C{. If the slice C{ 
is normalized, x remains constant: x(C J v ) = x ^03 f° r all normalizations 03. 

Therefore x leads to a strongly stable slice-sorting algorithm. The converse is also 
true if d = 2 , that is, if the slice is one-dimensional, then if 

h(cl) = h(m (ci)) 

for all normalizations 03 then h can only depend on the number of allocated {true) 
values in the slice since it fully characterizes the slice up to normalization. For the 
general case (d > 2), the converse is not true since the number of allocated values 
is not enough to characterize the slices up to normalization. For example, one could 
count how many sub-slices along a chosen second attribute have no allocated value. 

A function g is symmetric if go 03 ~ g for all normalizations 03. The following 
proposition shows that up to a monotone function, strongly stable slice-sorting al- 
gorithms are characterized by symmetric functions . 

Proposition 12 A slice-sorting algorithm based on a function g is strongly stable 
if and only if for any normalization 03, there is a monotone function h : IR — > R such 
that 

g(®(ci))=h(g(CJ)) (3) 

for all attribute values v = 1 , ... ,n of all attributes j = 1 , . . . , d. In other words, it 
is strongly stable if and only if g is symmetric. 

PROOF. By Proposition 10, Equation 3 is sufficient for strong stability. On the 
other hand, suppose that the slice-sorting algorithm is strongly stable and that 
there does not exist a strictly monotone function h satisfying Equation 3, then 
by Proposition 9, there must be a sorted sequence g{C J Vl ),g{C J V2 ),g{C J Vi ) such that 
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Table 2 

Examples of 2-d data cubes and their probability distributions. 



Data Cube 


Joint Prob. Dist. 


Joint Independent Prob. Dist. 
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g (w (ch^ '8 (® (^2)) (® (^ 3 )) * S n0t sortec ^ because this last statement 
contradicts strong stability, we have that Equation 3 is necessary. 



□ 



Lemma 13 A slice-sorting algorithm based on a function g is strongly stable if 
g = ho xfor some function h. For 2-d cubes, the condition is necessary. 

In the above lemma, whenever h is strictly monotone, then g ~ x and we call this 
class of slice-sorting algorithms Frequency Sort [9]. We will show that we can 
estimate a priori the efficiency of this class (see Theorem 18). 

It is useful to consider a data cube as a probability distribution in the following 
sense: given a data cube C, let the joint probability distribution *P over the same n d 
set of indices be 

'i/#cifq 1; ...,„^o 

otherwise 



The underlying probabilistic model is that allocated cells are uniformly likely to be 
picked whereas unallocated cells are never picked. Given an attribute j e {1 , . . . , d}, 
consider the number of allocated slices in slice C{, #Ci, for v G {1, . . . ,«}: we can 

define a probability distribution cp- 7 along attribute j as (pi, = From these (p J for 
all j e { 1 , . . . , d}, we can define the joint independent probability distribution $ as 



®h,-,id = IT/=i <P/., or in other words 4> = q>° ( 
Table 2. 



ty d 1 . Examples are given in 



Given a joint probability distribution *P and the number of allocated cells #C, we 
can build an allocation cube A by computing *P x #C. Unlike a data cube, an allo- 
cation cube stores values between and 1 indicating how likely it is that the cell 



17 



be allocated. If we start from a data cube C and compute its joint probability dis- 
tribution and from it, its allocation cube, we get a cube containing only O's and l's 
depending on whether or not the given cell is allocated (1 if allocated, otherwise) 
and we say we have the strict allocation cube of the data cube C. For an alloca- 
tion cube A, we define #A as the sum of all cells. We define the normalization of 
an allocation cube in the obvious way. The more interesting case arises when we 
consider the joint independent probability distribution: its allocation cube contains 
O's and l's but also intermediate values. Given an arbitrary allocation cube A and 
another allocation cube B, A is compatible with B if any non-zero cell in B has 
a value greater than the corresponding cell in A and if all non-zero cells in B are 
non-zero in A. We say that A is strongly compatible with B if, in addition to being 
compatible with B, all non-zero cells in A are non-zero in B Given an allocation 
cube A compatible with B, we can define the strongly compatible allocation cube 
A B as 



and we denote the remainder by Age = A — Ag. The following result is immediate 
from the definitions. 

Lemma 14 Given a data cube C and its joint independent probability distribution 
4>, let A be the allocation cube of<t>, then we have A is compatible with C. Unless 
A is also the strict allocation cube ofC, A is not strongly compatible with C. 

We can compute H(A), the HOLAP cost of an allocation cube A, by looking at each 
block. The cost of storing a block densely is still M = mi x . . . x whereas the cost 
of storing it sparsely is (d/2+l)D where D is the sum of the 0-to-l values stored 
in the corresponding block. As before, a block is stored densely when D > jjw+\\ ■ 
When B is the strict allocation cube of a cube C, then H(C) = H(B) immediately. If 
#A = #5 and A is compatible with B, then H(A) > H(B) since the number of dense 
blocks can only be less. Similarly, since A is strongly compatible with B, A has the 
set of allocated cells as B but with lesser values. Hence H(A) < H(B). 

Lemma 15 Given a data cube C and its strict allocation cube B, for all allocation 
cubes A compatible with B such that #A = #5, we have H(A) > H(B). On the 
other hand, if A is strongly compatible with B but not necessarily #A = #5, then 



A corollary of Lemma 15 is that the joint independent probability distribution gives 
a bound on the HOLAP cost of a data cube. 

Corollary 16 The allocation cube A of the joint independent probability distribu- 
tion 4> of a data cube C satisfies H(A)> H(C). 

Given a data cube C, consider a normalization 03 such that //(03(C)) is minimal and 




H(A) < H(B). 
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fs E Sc,x- Since H(fs(C)) < H(fs(A)) by Corollary 16 and //(03(C)) > #C by our 
cost model, then 

H(fs(C)) -//(03(C)) < H(fs(A)) -#C. 

In turn, H(fs(A)) may be estimated using only the attribute- wise frequency distribu- 
tions and thus we may have a fast estimate of H(fs(C)) —//(03(C)). Also, because 
joint independent probability distributions are separable, Frequency Sort is optimal 
over them. 

Proposition 17 Consider a data cube C and the allocation cube A of its joint in- 
dependent probability distribution. A Frequency Sort normalization fs G Sc,x is op- 
timal over joint independent probability distributions ( H(fs(A)) is minimal ). 

PROOF. In what follows, we consider only allocation cubes from independent 
probability distributions and proceed by induction. Let D be the sum of cells in 
a block and let Fa(x) = #(D > x) and /a(x) = #(D = x) denote, respectively, the 
number of blocks where the count is greater than (or equal to) x for allocation cube 
A. 

Frequency Sort is clearly optimal over any one-dimensional cube A in the sense that 
in minimizes the HOLAP cost. In fact, Frequency Sort maximizes Fa(x), which is 
a stronger condition (Ff s u) (x) > Fa(x)). 

Consider two allocation cubes A\ and A2 and their product A\ ®A2. Suppose that 
Frequency Sort is an optimal normalization for both A 1 and A2. Then the following 
argument shows that it must be so for A\ <g> A2. Block-wise, the sum of the cells in 

S\ S\ S\ S\ S\ 

A\ ® A2, is given by D = D\ x D2 where D\ and D2 are respectively the sum of 
cells in Ai and A2 for the corresponding blocks. 

We have that 

F Ai ®a 2 (x) = Y,fA ] (y)F Al (x/y) = Y,F A] (x/y)f A2 (y) 

y y 

and/s(Ai <g> A2) =fs(A\) ®fs(A2). By the induction hypothesis, Ff s f Al \(x) > Fa 1 (x) 
and so "L y F Ax {x/y)f Al {y) < L y Ffi{A 1 )(x/y)fA 2 (y). But we can also repeat the argu- 
ment by symmetry 

Y, F (f s ( A Mx/y)fA 2 (y) = Y,ffs(A 1 )(y) F A 2 (x/y) < Y.fHAi)(y) F fs{A 2 ){x/y) 

y y y 

and so Fa x ®a 2 ( x ) ^ F fs{Ai®A 2 ) ( x ) ■ The result then follows by induction. □ 

There is an even simpler way to estimate H(fs(C)) —//(03(C)) and thus decide 
whether Frequency Sorting is sufficient as Theorem 18 shows (see Table 3 for ex- 
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amples). It should be noted that we give an estimate valid independently of the 
dimensions of the blocks; thus, it is necessarily suboptimal. 

Theorem 18 Given a data cube C, let 05 be an optimal normalization andfs be a 
Frequency Sort normalization, then 

H(fs(C))-H(®(C))<(^ + l)(l-<l>-B)#C 

where B is the strict allocation cube ofC and 4> is the joint independent probability 
distribution. The symbol ■ denotes the scalar product defined in the usual way. 



PROOF. Let A be the allocation cube of the joint independent probability distri- 
bution. We use the fact that 

H(fs(C))-H(®(C)) < H(fs(A)) -//(03(C)). 

We have that fs is an optimal normalization over joint independent probability dis- 
tribution by Proposition 17 so that H(fs(A)) < //(05(A) ). Also //(05 (C) ) = //(G5(£) ) 
by definition so that 



H(fs(C)) -//(05(C)) < //(05(A)) -7/(05(5)) 

<//(G5(A B )) +H(U5(A B c)) -H(tH{B)) 
<//(G5(A se )) 



since H(tB(A B )) - //(05(A)) < by Lemma 15 
Finally, we have that 

H(m(Avc)\ < I _ 



H{US(A B c))< ( ^ + 1 )#A B c 



and#A sc = (1 -<&B)#C. □ 



This theorem says that 4> ■ B gives a rough measure of how well we can expect 
Frequency Sort to perform over all block dimensions: when $ • B is very close 
to 1, we need not use anything but Frequency Sort whereas when it gets close to 
0, we can expect Frequency Sort to be less efficient. We call this coefficient the 
Independence Sum. 

Hence, if the ROLAP storage cost is denoted by rolap, the optimally normalized 
block-coded cost by optimal, and the Independence Sum by IS, we have the rela- 
tionship 

rolap > optimal + (1 —IS)rolap >fs> optimal 
where fs is the block-coded cost using Frequency Sort as a normalization algorithm. 
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Table 3 

Given data cubes, we give lowest possible HOLAP cost H(US(C)) using 2x2 blocks, and 
an example of a Frequency Sort HOLAP cost H(fs(C) ) plus the independence product <1> B 
and the bound from theorem 18 for the lack of optimality of Frequency Sort. 
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6 Heuristics 



Since many practical cases appear intractable, we must resort to heuristics when the 
Independence Sum is small. We have experimented with several different heuris- 
tics, and we can categorize possible heuristics as block-oblivious versus block- 
aware, dimension-at-a-time or holistic, orthogonal or not. 

Block-aware heuristics use information about the shape and positioning of blocks. 
In contrast, Frequency Sort (FS) is an example of a block-oblivious heuristic: it 
makes no use of block information (see Fig. 3). Overall, block-aware heuristics 
should be able to obtain better performance when the block size is known, but 
may obtain poor performance when the block size used does not match the block 
size assumed during normalization. The block-oblivious heuristics should be more 
robust. 

All our heuristics reorder one dimension at a time, as opposed to a "holistic" ap- 
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input a cube C 

for all dimensions i do 

for all attribute values v do 

Count the number of allocated cells in corresponding slice (value of #C l v ) 

end for 

sort the attribute values v according to #C' V 
end for 

Fig. 3. Frequency Sort (FS) Normalization Algorithm 

proach when several dimensions are simultaneously reordered. In some heuristics, 
the permutation chosen for one dimension does not affect which permutation is 
chosen for another dimension. Such heuristics are orthogonal, and all the strongly 
stable slice-sorting algorithms in Section 5 are examples. Orthogonal heuristics 
can safely process dimensions one at a time, and in any order. With non-orthogonal 
heuristics that process one dimension at a time, we typically process all dimensions 
once, and repeat until some stopping condition is met. 



6. 1 Iterated Matching heuristic 



We have already shown that the weighted-matching algorithm can produce an op- 
timal normalization for blocks of size 2 (see Section 4.1.1). The Iterated Matching 
(IM) heuristic processes each dimension independently, behaving each time as if 
the blocks consisted of two cells aligned with the current dimension (see Fig. 4). 
Since it tries to match slices two-by-two so as to align many allocated cells in 
blocks of size 2, it should perform well over 2-regular blocks. It processes each 
dimension exactly once because it is orthogonal. 



This algorithm is better explained using an example. Applying this algorithm along 
the rows of the cube in Fig. 2 (see page 12) amounts to building the graph in the 
same figure and solving the weighted-matching problem over this graph. The cube 
would then be normalized to 

""l 1 - 

- 1 - 

1 

1 - 1 



We would then repeat on the columns (over all dimensions). A small example, 
1—1 1 1 

, demonstrates this approach is suboptimal, since the normalization 

1 

shown is optimal for 2 x 1 and 1x2 blocks but not optimal for 2 x 2 blocks. 
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input a cube C 

for all dimensions i do 

for all attribute values vi do 
for all attribute values V2 do 

w VljV2 <— storage cost of slices Q and CJ, using 

blocks of shape 1 x . . . x 1 x 2 x 1 x . . x 1 

i— 1 rf— (' 

end for 
end for 

form graph G with attribute values v as nodes and edge weights w 
solve the weighted-matching problem over G 

order the attribute values so that matched values are listed consecutively 
end for 

Fig. 4. Iterated Matching (IM) Normalization Algorithm 
6.2 One-Dense-Chunk Heuristic: iterated Greedy Sort ( GS) 



Earlier work [9] discusses data-cube normalization under a different HOLAP model, 
where only one block may be stored densely, but the block's size is chosen adap- 
tively. Despite model differences, normalizations that cluster data into a single large 
chunk intuitively should be useful with our current model. We adapted the most 
successful heuristic identified in the earlier work and called the result GS for iter- 
ated Greedy Sort (see Fig. 5). It can be viewed as a variant of Frequency Sort that 
ignores portions of the cube that appear too sparse. 

This algorithm's details are shown in Fig. 5 and sketched briefly next. Parameter 
pbreak-even can be set to the break-even density for HOLAP storage (pbreak-even = 
aJ+T = d/2+i ) ( see sec, i on 2)- The algorithm partitions every dimension's values 
into "dense" and "sparse" values, based on the current partitioning of all other 
dimensions' values. It proceeds in several phases, where each phase cycles once 
through the dimensions, improving the partitioning choices for that dimension. The 
choices are made greedily within a given phase, although they may be revised in a 
later phase. The algorithm often converges well before 20 phases. 

Figure 6 shows GS working over a two-dimensional example with pbreak-even = 
jw+l = \- The goal of GS is to mark a certain number of rows and columns as 
dense: we would then group these cells together in the hope of increasing the num- 
ber of dense blocks. Set A,- contains all "dense" attribute values for dimension /. 
Initially, A, contains all attribute values for all dimensions i. The initial figure is not 
shown but would be similar to the upper left figure, except that all allocated cells 
would be marked as dense (dark square). In the upper-left figure, we present the 
result after the rows (dimension i = 1) have been processed for the first time. Rows 
other than 1, 7 and 8 were insufficiently dense and hence removed from Ai: all al- 
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input a cube C, break-even density pbreak-even = d/2+1 
for all dimensions i do 

{A,- records attribute values classified as dense (initially, all)} 

initialize A,- to contain each attribute value v 
end for 

for 20 repetitions do 
for all dimensions i do 

for all attribute values v do 

{current A values mark off a subset of the slice as "dense"} 
p v <— density of C' v within Ai x A2 x . . . x A,-_i x A,-+i x . . . 
if p v < pbreak-even and v e A, then 

remove v from A, 
else if p v > pbreak-even and v ^ A,- then 

add v to Ai 
end if 
end for 

if A,- is empty then 

add v to A,, for an attribute v maximizing p v 
end if 
end for 
end for 

Re-normalize C so that each dimension is sorted by its final p values 

Fig. 5. Greedy Sort (GS) Normalization Algorithm 

located cells outside these rows have been marked "sparse" (light square). Then the 
columns (dimension i = 2) are processed for the first time, considering only cells 
on rows 1 , 7 and 8, and the result is shown in the upper right. Columns 0, 1,3,5 and 
6 are insufficiently dense and removed from A2, so a few more allocated cells were 
marked as sparse (light square). For instance, the density for column is i because 
we are considering only rows 1, 7 and 8. GS then re-examines the rows (using the 
new A2 = {2,4,7,8,9}) and reclassifies rows 4 and 5 as dense, thereby updating 
Ai = {1,4, 5, 7, 8}. Then, when the columns are re-examined, we find that the den- 
sity of column has become | and reclassify it as dense (A2 = {0,2,4,7,8,9}). 
A few more iterations would be required before this example converges. Then we 
would sort rows and columns by decreasing density in the hope that allocated cells 
would be clustered near cell (0,0). (If rows 4, 5 and 8 continue to be 100% dense, 
the normalization would put them first.) 

6.3 Summary of heuristics 

Recall that all our heuristics are of the type "1-dimension-at-a-time", in that they 
normalize one dimension at a time. Greedy Sort (GS) is not orthogonal whereas 
Iterated Matching (IM) and Frequency Sort (FS) are: indeed GS revisits the dimen- 
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Fig. 6. GS Example. Top left: after rows processed once. Top right: after columns processed 
once. Bottom left: after rows processed again. Bottom right: after columns processed again. 



sions several times for different results. FS and GS are block-oblivious whereas IM 
assumes 2-regular blocks. The following table is a summary: 



Heuristic 


block-oblivious/block-aware 


orthogonal 


FS 


block-oblivious 


true 


GS 


block-oblivious 


false 


IM 


block-aware 


true 
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7 Experimental Results 



In describing the experiments, we discuss the data sets used, the heuristics tested, 
and the results observed. 

7. 1 Data Sets 

Recalling that E[C) measures the cost per allocated cell, we define the kernel 
Km\,...,m d as the set of all data cubes C of given dimensions such that E(C) is min- 
imal (E(C) = 1) for some fixed block dimensions mi, . . . ,m^. In other words, it is 
the set of all data cubes C where all blocks have density 1 or 0. 

Heuristics were tested on a variety of data cubes. Several synthetic 12 x 12 x 12 x 
12 data sets were used, and 100 random data cubes of each variety were taken. 

• 2 2 re f ers to choosing a cube C uniformly from K2,2,2,2 and choosing n uni- 
formly from the set of all normalizations. Cube 71 (C) provides the test data; a 
best-possible normalization will compress 7t(C) by a ratio of max(p, |), where p 
is the density of 7t(C). (The expected value of p is 50%.) 

• K^ 2 2 2 is similar, except that the random selection from K2 ; 2,2,2 is biased towards 
sparse cubes. (Each of the 256 blocks is independently chosen to be full with 
probability 10% and empty with probability 90%.) The expected density of such 
cubes is 10%, and thus the entire cube will likely be stored sparsely. The best 
compression for such a cube is to j of its original cost. 

• K 2 2 2 2 + ^ a(lt ^ s n °i se - F° r every index, there is a 3% chance that its status (al- 
located or not) will be inverted. Due to the noise, the cube usually cannot be 
normalized to a kernel cube, and hence the best possible compression is proba- 
bly closer to | + 3%. 

• K 4^4 4 4+^ is similar, except we choose from 1^4,4,4, not K2,2,2,2- 

Besides synthetic data sets, we have experimented with several data sets used pre- 
viously [21]: Census (50 6-d projections of an 18-d data set) and Forest (50 
3-d projections of an 1 1-d data set) from the KDD repository [22], and Weather 
(50 5-d projections of an 18-d data set) [230 These data sets were obtained in 
relational form, as a sequence (t) of tuples and their initial normalizations can be 
summarized as "first seen, first when normalized," which is arguably the normal- 
ization that minimizes data-cube implementation effort. More precisely, let n be 
the normal relational projection operator; e.g., 

% 2 (((a,b),(c,d),(e,f))) = (b,d,f). 

2 Projections were selected at random but, to keep test runs from taking too long, cubes 
were required to be smaller than about 100MB. 
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Table 4 

Performance of heuristics. Compression ratios are in percent and are averages. Each num- 
ber represents 100 test runs for the synthetic data sets and 50 test runs for the others. Each 
experiment's outcome was the ratio of the heuristic storage cost to the default normaliza- 
tion's storage cost. Smaller is better. 



Heuristic 


Synthetic Kernel-Based Data Sets 


"Real- World" Data Sets 




^2,2,2,2 


K Sp 

^2,2,2,2 


sp N 

K 2,2,2,2 +N 


sp , T 


Census 


Forest 


Weather 


FS 


61.2 


56.1 


85.9 


70.2 


78.8 


94.5 


88.6 


GS 


61.2 


87.4 


86.8 


72.1 


79.3 


94.2 


89.5 


IM 


51.5 


33.7 


49.4 


97.5 


78.2 


86.2 


85.4 


Best result (estimated) 


40 


33 


36 


36 









Also let the rank r(v, (t)) of a value v in a sequence (t) be the number of distinct 
values that precede the first occurrence of v in (t). The initial normalization for a 
data set (t) permutes dimension i by y, where y^ 1 (v) = r(%i( (t) ) ) . If the tuples were 
originally presented in a random order, commonly occurring values can be expected 
to be mapped to small indices: in that sense, the initial normalization resembles an 
imperfect Frequency Sort. This initial normalization has been called "Order 7" in 
earlier work [9] . 



7.2 Results 

The heuristics selected for testing were Frequency Sort (FS), Iterated Greedy Sort 
(GS), and Iterated Matching (IM). Except for the "k^ 44 +N" data sets, where 4- 
regular blocks were used, blocks were 2-regular. IM implicitly assumes 2-regular 
blocks. Results are shown in Table 4. 

Looking at the results in Table 4 for synthetic data sets, we see that GS was never 
better than FS; this is perhaps not surprising, because the main difference between 
FS and GS is that the latter does additional work to ensure allocated cells are within 
a single hyperrectangle and that cells outside this hyperrectangle are discounted. 

sp sp 

Comparing the K 2 2 2 2 and K 2 2 2 2 +N columns, it is apparent that noise hurt all 

heuristics, particularly the slice-sorting ones (FS and GS). However, FS and GS 

sp sp 
performed better on larger blocks (k 4 4 4 4 +N) than on smaller ones (k 2 2 2 2 + N) 

whereas IM did worse on larger blocks. We explain this improved performance for 

slice-sorting normalizations (FS and GS) as follows: #C' V is a multiple of 4 3 under 

K4 5 4 5 4 5 4 but a multiple of 2 3 under K2,2,2,2- Thus, K2,2,2,2 is more susceptible to noise 

than K4,4,4,4 under FS because the values #C l v are less separated. IM did worse on 

larger blocks because it was designed for 2-regular blocks. 

Table 4 also contains results for "real-world" data, and the relative performance 
of the various heuristics depended heavily on the nature of the data set used. For 
instance, Forest contains many measurements of physical characteristics of geo- 
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Fig. 7. Solution-size ratios of FS and IM as a function of Independence Sum. When the 
ratio is above 1.0, FS is suboptimal; when it is less than 1.0, IM is suboptimal. We see that 
as the Independence Sum approached 1.0, FS matched IM's performance. 

graphic areas, and significant correlation between characteristics penalized FS. 



7.2.1 U tility of the Independence Sum 

Despite the differences between data sets, the Independence Sum (from Section 5) 

size usins FS 

seems to be useful. In Figure 7 we plot the ratio -■ ■ & T1V . against the Indepen- 

° v size using IM ° v 

dence Sum. When the Independence Sum exceeded 0.72, the ratio was always near 
1 (within 5%); thus, there is no need to use the more computationally expensive 
IM heuristic. Weather had few cubes with Independence Sum over 0.6, but these 
had ratios near 1.0. For Census, having an Independence Sum over 0.6 seemed 
to guarantee good relative performance for FS. On Forest, however, FS showed 
poorer performance until the Independence Sum became larger 0.72). 



7.2.2 Density and Compressibility 

The results of Table 4 are averages over cubes of different densities. Intuitively, for 
very sparse cubes (density near 0) or for very dense cubes (density near 100%), we 
would expect attribute- value reordering to have a small effect on compressibility: 
if all blocks are either all dense or all sparse, then attribute reordering does not 
affect storage efficiency. We take the source data from Table 4 regarding Iterated 
Matching (IM) and we plot the compression ratios versus the density of the cubes 
(see Fig. 8). Two of three data sets showed some compression-ratio improvements 
when the density is increased, but the results are not conclusive. An extensive study 
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Fig. 8. Compression ratios achieved with IM versus density for 50 test runs on three data 
sets. The bottom plot shows linear regression on a logarithmic scale: both CENSUS and 
Weather showed a tendency to better compression with higher density. 



of a related problem is described elsewhere [9]. 



7.2.3 Comparison with Pure ROLAP Coding 

To place the efficiency gains from normalization into context, we calculated (for 
each of the 50 Census cubes) Cdefauit, the HOLAP storage cost using 2-regular 
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blocks and the default normalization. We also calculated crolap, the ROLAP cost, 
for each cube. The average of the 50 ratios Cdefault was 0.69 with a standard devi- 
ation of 0.14. In other words, block-coding was 31% more efficient than ROLAP. 
On the other hand, we have shown that normalization brought gains of about 19% 
over the default normalization and the storage ratio itself was brought from 0.69 to 
0.56 in going from simple block coding to block coding together with optimized 
normalization. Forest and Weather were similar, and their respective average 
ratios Cdefault were 0.69 and 0.81. Their respective normalization gains were about 

crolap r ° 

14% and 12%, resulting in overall storage ratios of about 0.60 and 0.71, respec- 
tively. 



8 Conclusion 

In this paper, we have given several theoretical results relating to cube normaliza- 
tion. Because even simple special cases of the problem are NP-hard, heuristics are 
needed. However, an optimal normalization can be computed when 1x2 blocks are 
used, and this forms the basis of the IM heuristic, which seemed most efficient in 
experiments. Nevertheless, a Frequency Sort algorithm is much faster, and another 
of the paper's theoretical conclusions was that this algorithm becomes increasingly 
optimal as the Independence Sum of the cube increases: if dimensions are nearly 
statistically independent, it is sufficient to sort the attribute values for each dimen- 
sion separately. Unfortunately, our theorem did not provide a very tight bound on 
suboptimality. Nevertheless, we determined experimentally that an Independence 
Sum greater than 0.72 always meant that Frequency Sort produced good results. 

As future work, we will seek tighter theoretical bounds and more effective heuris- 
tics for the cases when the Independence Sum is small. We are implementing the 
proposed architecture by combining an embedded relational database with a C++ 
layer. We will verify our claim that a more efficient normalization leads to faster 
queries. 
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